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Abstract 

We present the exact 0(a) correction to the process e + e~ — > ff + 7, / 7^ e, 
for ISR©FSR at and beyond LEP2 energies. We give explicit formulas for the 
completely differential cross section. As an important application, we compute the 
size of the respective sub-leading corrections of O(aL) to the // cross section, where 
L is the respective big logarithm in the renormalization group sense so that it is 
identifiable as L = In |s|/m e 2 when s is the squared e + e~ cms energy. Comparisons 
are made with the available literature. We show explicitly that our results have 
the correct infrared limit, as a cross-check. Some comments are made about the 
implementation of our results in the framework of the Monte Carlo event generator 
OC MC. 
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Currently, the final LEP2 data analysis is in its beginning stages, and the desired 
total precision tags on the important LEP2 physics processes e + e~ — > //, / ^ e, are 
already called out in the LEP2 MC Workshop in Ref. |]]. It has been demonstrated in 
Ref. that the Monte Carlo (MC) event generator program /C/C f2|, hereafter referred 
to as /C/C MC , and the semi-analytical program ZFITTER || realize these precisions 
(.2 — 1%) in most channels for inclusive cross sections and that for the fully differential 
distributions, the /C/C MC again meets most of the requirements for the LEP2 final data 
analysis. In this paper, we present exact results on the 0(a) correction to the single 
hard bremsstrahlung processes e + e~ — ► // + 7, / ^ e. This correction is an important 
contribution to the differential distributions as they are realized in the /C/C MC which 
allows the very demanding precisions just cited to be achieved. 

Specifically, the exact results for the 0{a) corrections to s-channel annihilation hard 
bremsstrahlung processes under study here were also considered in Refs. We differ 

from these results as follows. Concerning Ref. 0, the entire result was given only for the 
case in which the photon angle variables are all integrated out; here, we give the fully 
differential results. With regard to Ref. 0], the completely differential results were given 
as well but the mass corrections were omitted. In our work, the masses of the electrons 
and positrons and the masses of the final state fermion and anti-fermion are exactly taken 
into account in contrast to the literature. Thus, by comparing with the two calculations 
in Refs. P,H as we do here, we get a measure of the size of the mass corrections as well 
as cross checks on both our differential and our integrated results. 

Our work is organized as follows. In Section [I], we set our notational conventions. 
In Section we present our exact amplitudes for the 0(a) virtual corrections to initial- 
state and final-state real radiation. In Section |3], we derive the differential cross-sections 
corresponding to these amplitudes in a form useful for comparisons. In Section ^, we 
compare these results with those in Refs. [||, [J while illustrating our results as they are 
used in the /C/C MC in Ref. ||. Section ^] contains our summary remarks. The Appendix 
contains technical details about the scalar integrals. 

1 Preliminaries 

In this section we set our notational conventions. We will use the conventions of Refs. 0, 
||,[7j] for our spinors. These conventions are based on the Kleiss-Stirling || Weyl spinors 
augmented as described in Refs. @,||[7j with the rules for controlling their complex phases, 
or equivalently, the three axes of the fermion rest frame in which the spin of that fermion is 
quantized. We sometimes refer to this fermion rest frame as the global positioning of spin 
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(GPS) frame and to the rule for determining it as the GPS rule. The resulting conventions 
for the fermion spinors are then called the GPS spinor conventions. See Refs. [@,|6|,[7| for 
more details. Let us now turn to the kinematics. 





v\,. y e > 



Figure 1: Representative initial state radiation graphs for e + e — > ff with one virtual and 
one real photon, for / ^ e. 

The process under discussion here is shown in Fig. [I], the one-loop virtual correction to 
the hard bremsstrahlung process e~e + — ► ff + j, for / 7^ e. We will treat both the initial 
state radiation (ISR) case and the final state radiation (FSR) case. We denote the four 
momenta and helicity of the e~, e + , /, and / as pj and Xj, j = 1, ...,4, respectively. We 
denote the charge of / by Qf in units of the positron charge e. The rest mass of fermion 
/ is denoted by nif. The photon momentum and helicity will be denoted by k and a. 
With our GPS conventions for spinors, we induce the following polarization vectors for 
photons: 



1.1 



y/2 u_ a (k)u a ((3) ' ffVbyy y/2 u_ a (k) U(7 (0 ' 
where the auxiliary 4-vector (3 is exploited here to simplify our expressions as needed. 
It satisfies j3 2 = 0. The second choice with u a (C), as defined in Ref. 0, is already an 
example of this exploitation - it often leads to simplifications in the resulting photon 
emission amplitudes. 

The calculations which we present have been done using the program FORM of Ref. . 
For the t-channel case, we presented similar results in Ref. [l(J in connection with the 
respective 0(a 2 )L corrections needed for the 0.061% (0.054%) total precision tag achieved 
in Ref. JTT| for the LEP1/SLC luminosity process in the Monte Carlo event generator 
BHLUMI4.04 in Ref. ||12|| . Just as in the latter case, here a considerable effort is needed 
to simplify our initial raw FORM output in order to make a practical application of the 
respective results in the context of a Monte Carlo environment such as the /C/C MC in 
Ref. 0. As in Ref. fTD"! , we only present the final simplified expressions in this paper for 
the sake of clarity. 

Our metric is that of Bjorken and Drell in Ref. [|13|] and we effect our gauge invariant 
calculation in the 't Hooft-Feynman gauge. With these preliminary remarks, we turn now 
in the next section to the calculation of our process of interest. 
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2 Exact Results on the Virtual Correction to e + e — ► 

ff + 7, f?e 

In this section we calculate the exact virtual correction to e + e~ — ► // + 7, f ^ e. We 
proceed in analogy with our results on the virtual correction for the t-channel dominated 
low angle Bhabha scattering process with a single hard bremsstrahlung in Ref. ||10|| . 

Specifically, we express the exact amplitude for one real and one virtual photon emitted 
from the electron lines in the process e + e _ — > f f + 7 using the GPS conventions [@,||,[7|]- 



In Ref. ||1 0|| , the corresponding t channel result was obtained for electron line emission. 
Here, from the latter result, we first obtain the respective initial-state s channel result 
by crossing the outgoing electron line with the incoming positron line, and replacing 
the respective final state by //, while adding also Z boson exchange. The results are 
translated into GPS conventions. Then, in subsection |2.2|, we provide the detailed form 



factors appearing in the initial state amplitudes. The corresponding final state amplitudes 



are presented in subsections |27| and 2A . 



2.1 ISR s-Channel Exact Result 

In this subsection, we define notation and set up the exact contribution for one real photon 
and one virtual photon emitted from the electron lines in the process e + e~ — > // + 7. The 
amplitude for real plus virtual photon emission from the initial state may be written 



M m» = $y M fm (/o + M + h) (21) 

107T Z 



where the real photon emission amplitude is A^i SR ^°' ) , and the factors /o,i,2 contain spinor 
dependence. They will be specified in the next section. 

too (n\ 

In GPS conventions, the amplitude Ai 1 for the initial state radiation of a single 
photon is given by 



V 



A 



eQ 



2kp 



+ ^v( P2 ,X 2 )f a (- 7 ! 2 + m+ / k)M 1 u( Pl ,X l ) , (2.2) 
2kp 2 



where 



Ml = ^^Q S y ( "(P3^3)W- B ^(P4,A 4 ; ) B ( } 

1 Z-~/ \ o/_|W2 1 ,T J u / ^ V ; 



s' - M% + iYbs'/Mb 
is the annihilation scattering spinor matrix. 



A,,j=±l 
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The form factors may be obtained from the corresponding t-channel result in Ref. [10 
for electron line emission. Specifically, the s channel result can be obtained by crossing 
the outgoing electron line with the incoming positron line, and replacing the final state 
by //• We also include the effects of Z exchange in the s channel. 



Our previous calculations of t-channel bremsstrahlung |T(| used the Chinese Magic 
conventions for the photon polarizations. The GPS version of the magic photon 
polarization vector is related to the Chinese Magic conventions [l4j by 



>GPS 



ae 



Chinese 



-cr 



(2.4) 



The purpose of this change is to recover the more standard convention of defining photon 
polarization in terms of incoming states. The choice of magic polarization vector affects 
the amplitude ( |2.1D only through the definition of A^i SR ^ ' ) . The remaining factors may 
thus be obtained directly by crossing from our previous t-channel results. 

The magic choice of auxiliary vector for initial state radiation is (3 = h 
the definition 



p 2 pi 

a Ai A2 



with 



h 



Qo <?1 <?2 

Mo A*i M2 







' Mi 


= M2 




H 


Mo 


= Mi 


?2 , 




k Mo 


= M2 



(2.5) 



"Mi 



Using the magic polarization vector in ( |2.1| ) and neglecting fermion masses gives 

,ISR(0) _ ,^ _„3^y fJ\ T 2s_ CT (p 3 ,p 4 ) 



iQ e ae G Xl ,\ 3 (s')I c 



s a (p 1 ,k)s a {p 2 ,k) 



(2.6) 



where the photon-Z propagator is 

GxM = 

and 



/ 



e,B f B 



B=-y,Z 



s' - M 2 B + iT B s' /M B ' 



(2.7) 



p 2 Pi 
a Xi X 2 



h 



0p 3 p 4 
a A 3 A 4 



(2.8) 



J = -V^X^sl lh 
We now turn to calculating the form factors and spinor factors. 

2.2 Initial State Form Factors 

It remains to describe the form factors and spinor factors needed to compute M. x ' . 
The spinor factors I\ 2 are given by 
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h = \ / 2\ 1 s_ Xl (p 1 ,k)s Xl (p 2 ,k) 

y S- Xl (P^Pl) s \i(PuP3) ~ S- Xl (p 4 ,p 2 )s Xl (p 2 ,p 3 ) ^ g * 

J2 = V2as- Xl (Pu k)s Xl (P2, k)s- X[i (p 4 , k)s Xs (p 3 , fe) ^ ^ 

S- CT (pi,p 2 )S- CT (p4,P3)^0 

where the spinor product is s x (p,q) = u^ x (p)u x (q). The factors Ii j2 are crossed versions 
of (Jj ± X 2 )/2J in Ref. 0. 

We will begin by writing the dominant term /q. Expressions can be found in Ref. || 
for all of the scalar integrals needed for the form factors, which were previously calculated 
using the FF package ||, which implements the methods of Ref. ||15|| . The integrals in 
Ref. || are not quite adequate, because of the possibility that < m 2 e / s. However, it was 
possible to analytically continue when necessary, and to reproduce the numerical results 
of the FF package. Thus, an expression for the form factors in terms of logarithms and 
dilogarithms is now available. Details on the s channel version of the scalar integrals used 
in Ref. Jill may be found in the Appendix. 

For a — Ai, using = 2p, • k/s, 

fo = 47rS YFS (s,m e ) + 2 (L - 1 - m) + ' 2 



1 — r 2 



+ / i n+ ? ) R(n,r 2 )+R(r 2 , ri ) 
(1 - r 2 )r 1 

+ i 3 + 71 XT— -rjlna-n-ra) 

I (1 -r 2 ){r 1 + r 2 ) J 

r2(2 + ri) V (1 - ri " r2) -z^ (2.11) 



(1 - n)(l - r 2 ) [ r 2 
with L = ln(s/TOg), the infrared YFS factor 

4^5^(3,771)= f41n^ + l) (ln^--l-Z7r)-ln 2 (^-)-l+^-+m(21n^--l) (2.12) 
and 

R(x,y) = ln 2 (l -x) + 21n(l -x) jln ^ +ivr 

+ 2Sp(x + y)-2Sp( T ^-) (2.13) 

= ln 2 (l -x) + 21n(l -x) jln (j - ^ + ivr 

+ 21npU lnf 1 - X -^ -2Sp^ * 



x + yj \ l — x J \x + y 

(1 — x — y)x 
(x + y)(l-x) 



2 Sp(x) + 2 Sp ( Z 1 \- V)X . ) . (2.14) 
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The second expression is preferred for calculating R(x,y)/x when x may be small. For 
a = —A, 7*1 and r 2 are interchanged in Q2.ll ). 



The coefficients of the spinor terms in fl2.1|) are, for a = A, 

(n-r 2 ) r 2 (l-n-r 2 ) 
71 2(l-r 1 )(l-r 2 )^n(l-r 2 )(ri+r 2 ) 1 1 2j 
ri + r 2 / 1 - n - r 2 1 \ n + r 2 



2ri [ ri(l — r 2 ) 2 ' J 4r 2 
1 — 7*i — r 2 ( r 1 + r 2 1 r 2 



and 



(1 — ri)(l— r 2 ) [2(1 — ri) 1 — ri r x 

x An 1 ~ ^ ~ r2 - ivr^ (2.15) 



, 9 2 - n - r 2 



2(1-7x1(1-7*2) 

l-r x -r 2 /2-r 2 V 

+ — ■ r ln(l-n-r 2 ) 

ri(ri + r 2 ) \1 - r 2 / 

2(1 - n - r ) ( i 

+ ^ + — — ln(l-n-r 2 

r*i + r 2 I n + r 2 



(1 _ r r ) (2+ri - r2) 

2rf(l - r 2 J 



+ ^ (i - r 2 ) - \ (i - ^ ] ^-^(r,, n) 

ri — r 2 {2-r 2 ^ r 2 - t*i 



(i_ ri )(i_ r2 ) ^ ri 2(1 -n) 



xlln- —-in). (2.16) 



r 2 

The coefficients /i )2 are s-channel versions of (JF X ± T^)!^ in Ref. [10]. For a = —A, 7*1 
and r 2 are interchanged in ( |2.15|) and ( |2.16| ). 

The leading log limit is obtained by finding which terms give rise to the leading powers 
of the 'big logarithm' L when the above expressions are integrated over r\ and r 2 . These 
come from collinear terms where 7*1 or r 2 go to zero. In the collinear limits, when averaged 
over the azimuthal angle, only the fo terms remain to order L 2 and L, i.e., to order NLL. 
Using the identities 



R(0,y) = 0, 

-R(x,y) = 2 

x v y 

R(x,y) = 21n(l -x) [hxy + iit) -ln 2 (l - x) + 2Sp(x) for y -> 0, (2.17) 



-R(x,y) = 2 ( 1 - - ) ln(l -y) - 2 my - 2ni for x -»• 



2/ 
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the NLL limit of the form factor /o is found to be 

jNll = 2(L- l-?7r) +21n(l-ri)(lnr 2 + i7r) +21n(l-r 2 )(lnri + i7r) 

- ln 2 (l - n) - ln 2 (l - r 2 ) + 31n(l - r x ) + 3 ln(l - r 2 ) 

+ 2 Sp(n) + 2 Sp(r 2 ) + 5 CT) _ Al + T^-^Aa (2.18) 

1 — T\ 1 — r 2 

without mass corrections. 

Mass corrections we have calculated primarily without any approximations, however, 
in the following we shall present them in the approximation m e << \/s. In particular, in 
this approximation, we checked by explicit calculation that the result which we obtain for 
the mass corrections in fact agrees with that implied by the prescription in Ref . ]TB| . This 
prescription is valid for the spin-averaged differential distribution in the limit m e << y/s, 
but since mass terms are located in the separate (helicity conservation violating) spin 
amplitudes, it is not difficult to "undo" the spin summation. The technique of Ref. fTE| 



was originally applied to tree level photon emissions. Following the Appendix B of Ref. [10 
we can apply it also to our case of emission of one virtual and one real photon. 

Taking advantage of the freedom, which we have for presenting mass terms in the 
m e << y/s approximation, the introduction of the mass correction leads to a replacement 
of /oby/ + /o me , where 

, (me) 2m e 2 fn t r 2 \ (l-ri)(l-r 2 ) 



s \r2 ' T\) 1 + (1 — T\ — r 2 ) 2 
x {/ -4vr J B YFS (s / ,m e )-2[L + ln(l-r 1 -r 2 )-l-m]} (2.19) 

with YFS infrared factor 

An B Y fs{s' , m e ) = Att B Y Fs{s,m e ) + ln(l - r 1 - r 2 ) (41n — 

m e 

- 2L + 1 + 2ni) - ln 2 (l - r\ - r 2 ) . (2.20) 
Mass corrections first appear at order NLL, and to this order, 

c(rn e ) NLL 2m e 2 ( T\ , T 2 \ (l-ri)(l-r 2 ) 



r\"i, e ) j\±j±j a i ■ x 



r 2 rij 1 + (1 — ri - r 2 



\2 



x |ln(l-ri-r 2 ) ( 2L - 1 - 2ni - 41n— J + ln 2 (l - n - r 2 ) 

k \ TTl e J 

+ 21n(l-ri)(lnr 2 + i7r)+21n(l-r 2 )(lnri+m)|. (2.21) 



Only the LL part of fo contributes to the mass correction, to order NLL. The result ( |2.19| ) 



gives the complete effect of the mass corrections for the ISR neglecting the terms that are 
suppressed by higher powers of m 2 e / s as usual. 
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2.3 Final State Radiation 



The amplitudes for final state radiation can be obtained by crossing the incoming electron 
with the outgoing /, and the incoming positron with the outgoing /. Thus, p\ <-» — p±, 
P2 —P3, Ai <-> — A 4 , and A 2 <-» — A 3 in the results of the previous sections. 

The final state radiation (FSR) amplitude can be written in analogy with the ISR 
result (O), 



M\ 



FSR(l) 



167T 2 ' 



^i SR(0) (/o + flh + f 2 h) 



(2.22) 



The form factors / 1 2 an d the spinor factors are final-state analogs of those in the 
previous section, and will be defined in the next subsection. 

The amplitude .A/fi SR( ' ^ for the final state radiation of a single photon can be obtained 
from the initial state amplitude M 1 ^^ by crossing. Crossing leads to spinors with 
negative energy. A consistent choice of branches gives 



u(-p, -A) = iv(p, A), v(-p, -A) = iu(p, A). 



Then we obtain 



FSR(O) 



e Q f _ 

7^M(P3,A 3 )Mi(y 4 -m + /k)£v(p4, A 4 ) 
eQ f 



2kp, 



u(p3, A3) £(y 3 + m + ^)Mi«(p 4 , A 4 



where 



Mi 



^(P2, A 2 )7^' B ^am(Pi,Ai 



B=7,Z 
A,/i=±l 



The magic polarization vector for final state radiation is h 
Q2.24 ) gives, in the massless limit, 

FSR(O) _„ 3 ^ 2s -AP3,P4.) 



P3 PA 
a A3 A4 



M[ 



%Q s ae i G. 



with propagator (277) and J given again by (|2.8|) . 



2.4 Final State Form Factors 



The spinor factors I\% appearing in are given by 

li = \/2A 4 sa 4 (p 4 , k)s_ Xi (P3,k) 

g A 4 (Pi, P4)g-A 4 (P2, + gA 4 (Pi, P3)g-A 4 (P3, jfr) 
S-<r(p4,P3)S- (T (p2,Pl)sA 4 (pi,P3)S-A4(P3,P2)i"o 

- \/2crSA4(P4, k)s_x 4 (P3, k)s X2 (pi, fc)s_ A2 (p 2 , fc) 

S_ <T (p4,P3)S- (T (pi,p 2 )/o 
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[2.23) 



(2.24) 

(2.25) 
Using this in 

(2.26) 



(2.27) 
(2.28) 



As before, let r, = 2pi ■ k/s. We can obtain the form factor f for a = — A 3 by substi- 
tuting 

n -»• -r 3 /(l -r 3 -r 4 ), r 2 -> -r 4 /(l - r 3 - r 4 ) (2.29) 

in (|2.11|) , and for a = +A3 by interchanging r 3 and r 4 in (|2.29|) . (Since A; 2 = 0, we have 
r i + r 2 = r 3 + r 4 .) Then, for cr = A 3 , 

m f . \ r 3 



1 — r 4 



/ = 47r J BYFs(s / ,m f ) + 2(L-21n^-l-i7r] 

V m e J 

(r 3 + r 4 )(l — r 3 — r 4 )— — 

+ v R{r3, U) + R{r 4 , r 3 ) 

r 4 (l - r 4 ) 

2r 3 (l - r 3 - r 4 ) \ 
+ ^5_/Jn_- 2 llnr 3 (2.30) 



1 — r 4 I 1 — r 3 



with 



-y -a; 



1 — x — 1 — £ — y 



= 21na;ln( — — — ) + 2 Sp(x) - 2 Sp(x + y). (2.31) 
Vl-x-y/ 

The expression ( j2.31| ) is obtained from ( [2.13 ) using the dilogarithm identity 



The imaginary parts in Q2.30T ) were obtained by assuming the «7r terms in ( p. 11 ) came 
from a small positive imaginary part on s or s' . For a = — A 3 , fo = /b( r 4; r 3) instead: r 3 
and r 4 are interchanged. 

The coefficients of the spinor terms in (|2.22|) are, for a = A 3 , 

- (r 3 -r 4 )(l-r 3 -r 4 ) r 3 (l - r 3 - r 4 ) 

2(1 - r 3 )(l - r 4 ) r 4 (l - r 4 )(r 3 + r 4 ) 
r 3 + r 4 Jl-r 3 -r 4 1 1- r 3 + r 4 — 

+ n (1 "?n" r4 \ + ~ + STpM ln (2 ' 33) 

(1 - r 3 )(l - r 4 ) [ r 4 2(1 - r 3 ) J 



-7 2 _ (1 - r 3 - r 4 )(2 - r 3 - r 4 ) 

2 2(l-r 3 )(l-r 4 ) 
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1 - r 3 - r 4 / r 3 



r 4 (r 3 + r 4 ) \ 1 - r 4 / 
2 f 1 — r 3 — r 4 



^ ln(l — r 3 — r 4 
ln(l — r 3 — r 4 ) 

r 3 + r 4 [ r 3 + r 4 

(1 -r 3 -r 4 )(2-r 3 -3r 4 ) — 

+ 2^r^) ^ (r3 ' r4) 

+ J (l - ^) ViEfo, r 4 ) - 1 (l - ^ ^(r,, r 3 ) 

r3 " r4 f2 ~ r3 + y 4 '" 3 - 2 llnr 3 . (2.34) 



(1 - r 3 )(l - r 4 ) [ r 4 2(1 - r 3 

For cr = — A 3 , r 3 and r 4 are interchanged in fl2.33j ) and ( |2.34|) . 

The NLL limit is obtained as in the initial state radiation case, except that now the 
collinear limits are when r 3 or r 4 become small. Only the form factor f survives to order 
NLL, and using the identities 



R{x,0) = 0, 
R(x,y) = - 

R(x,y) = -21na;ln(l-y)-2Sp(y) for x ^ 0, (2.35) 



1— 2 2 

-R(x,y) = lnxH — ln(l — x) for y — > 

y 1 — x x 



we find 



/^ iL = 47rB Y Fs(s',m / ) + 2(L-l-*7r) 

- 2 lnr 3 ln(l — r 4 ) — 2 lnr 4 ln(l — r 3 ) 

- ln(l - r 3 ) - ln(l - r 4 ) - 2 Sp(r 3 ) - 2 Sp(r 4 ) 

- ^ CT ,A 3 ^3 - 5(7 ,-A 3 »"4 (2.36) 

without mass corrections. 

Spin-averaged mass corrections can be obtained from the initial state case (|2.19|) by 
crossing. The result is that / — > f Q + f Q s , where 

(m/) 2mJ/r 3 r 4 \ (l-r 3 )(l-r 4 ) 



fo = — f + 



where 



Ja ' r 3 J 1 + (1 - r 3 — r 4 ) 2 
x |/ -47r J B Y FS/(s)-2(L-21n^-l-m)| (2.37) 



47r J B Y FSf(s) = 47r J B Y Fs(s',m / ) - ln(l - r 3 - r 4 )(41n — 

- 2L + 1 + 2ixi) + ln 2 (l - r 3 - r 4 ) . (2.38) 
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Again, mass corrections first appear at order NLL, and to this order, 
-p(mf) nll _ 2m| / r 3 r 4 ^ (l-r 3 )(l-r 4 ) 



s \r 4 r-^j 1 + (1 — r 3 — r 4 ) 2 

x / ln(l - r 3 - r 4 ) ( 4 In — - 2L + 4 In ^ + 1 + 2m } 
I \ mf m e J 

- In 2 (l-r3-r 4 )+21n(l-r 3 )lnr 4 + 21n(l-r 4 )lnr 3 |. (2.39) 

The result (|2.37| ) gives the complete effect of FSR mass corrections neglecting terms 
suppressed by higher powers of — as usual. 



3 Differential Cross Section 

This section translates our amplitudes into differential cross sections, and sets up com- 
parisons with other related results. The initial state differential cross section for emitting 
one real photon may be written 

j isr(o) . 9 

d(T i = 1 V U/f ISR(0) " (3 1) 

d 2 Qdndr 2 2(4tt)V ^ I 1 ' 1 ' 

where the summed, squared real photon amplitude leads to 

Mf m 2 = [(tl + tl)(F - F 1 ) + {u\ + ul)(F + Fl )] , (3.2) 

where the invariants Ui, ti may be written in YFS3-style ||[17| effective angle notation: 



U = iPi ~ Pi+2) 2 = --/3ff3is(l - Ti) [2 - cos(6 l ii) - cos(6»2i)] 
Ui =(pi-p j+2 ) 2 = --PffiisQ. -r i )[2 + cos(6 l H) + cos(6> 2i ) ] (3.3) 



4 

with (ij) = (1,2) or (2,1), and 



4m 2 f / 4m 2 s' 

We will be setting these mass factors to unity in the following, and adding mass corrections 
at the end via ([TTJ) or (|Z37|). 



The coefficients Fi are defined in terms of the standard vector and axial vector fermion 
couplings V and A, and X = V 2 + A 2 , Y = 2VA by 

F = X e X fX 2 + 2Q e Q f V e V f Xi + Q 2 e Q 2 f , 

Fi = Y e Y fX 2 + 2Q e Q f A e A fXl , (3.5) 



12 



where 



Xi = s'{s'-Ml)[{s'-M 2 z f + {s'V z /M z f]-\ 

X 2 = A(s' -M 2 z f + {s'T z /M z f}-\ (3.6) 

The initial state differential cross section for real plus virtual photon emission may be 
expressed as 

j ISR(l) , 

= -r^— T Re UM?Wy M *W)-\ . ( 3. 7) 
d 2 ndr 1 dr 2 (47r)V ^ I J 

Xj,CT 

If is convenient to rewrite (12.11) as 



M mm = gy v m isr(o) (38) 

Ion 2 

in terms of a virtual correction factor 

« = fo + fih + hh- (3.9) 

The differential cross section ( 3.7|) can then be written in terms of a spin-averaged virtual 
correction factor (w) times the cross section for pure real initial state radiation: 

daf R{1) = (v) daf R{0 \ (3.10) 

where 

£^ 3 ^! SR(0) I 2 



(v) 



E Al Kl^o + + Ai(F 2 + F 3 )] + v^[F -F 1 + Ai(F 2 - F 3 )}} 
(Fo + F^uj + uD + iFo-F^tl + tl) 



(3-11) 



with 



and 



= A,A -A + «^A,A,A> 

V- = iK-A,-A + *H-A,A, (3.12) 



F 2 = F e X /X 2 + 2g e g / A e V/Xi, (3.13) 
F 3 = X e Y fX 2+2QeQfV e A f xt. 
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In the NLL approximation, where the f\ and f 2 terms in (|2.1| ) may be neglected, we 
may use the relation (v ) = (fo)- Then Q3.ll ) simplifies to 

/ f \ _ o (/o(ri,r 2 )^+/o(r 2 ,ri)^)(F +Fi)+(/ (ri r 2 )f|+/ (r 2 ,ri)ff)(F -Fi) 
WO/ z ( u 2 ,'2 s? F T p T ; .2 T 1 3 w Fn _ pT s 



K +« 2 ) (Fo +Fi ) + (t? +^ ) (-Fo - fi ) 
Fq[/ (''1i''2)(1+cos 2 6>n)+/ (r 2 ,ri)(l+cos 2 6 22 )]+2Fi [f (ri ,r 2 ) cos 6>n+/ (r 2 ,ri) cos 6> 22 ] 



(3.14) 



F [(l-ri) 2 (l+cos 2 6»n)+(l-r 2 ) 2 (l+cos 2 6> 22 )]+2Fi [(1-n) 2 cos 0n+(l-r 2 ) 2 cos 6» 22 ] 

with 

/ (r l ,r,) = 2(l-r,) 2 /o(r i ,r,). (3.15) 

If we further drop the dependence on 6^ in ( |3.3|) , letting 9ij — > 9 for fixed 9, then we get 
a simpler approximation. The angle dependence can be factored out of the cross-section 
leading to 

At! SR(1) 1 /Q e e 2 1 



dr 1 dr 2 2 V 2tt 



ao(f )H ( ri ,r 2 ). (3.16) 



with the definition 



Ho(n, r 2 ) = - [(1 - n f + (1 - r 2 ) 2 ] , (3.17) 

lr\r 2 



the total Born cross section 



cr = ^ e ( / /rf 2 -F (l + cos 2 0)+FiCOS0), (3.18) 



,, x I-l - ' 2; 70V 1, ' 2) T ^ - ' U JOV 2, ' 1; /QirA 
Uo) - 7^ _ yji n _„.V2 • ^- iy J 



2(4vr) 2 s'7 V 2 " 
and approximate spin-averaged form factor 

(1 - r 2 ) 2 f Q {r l ,r 2 ) + (1 - ri) 2 / (r 2 , r x ) 

;i-n) 2 + (l-r 2 ) 

It can be shown that this approximation is valid to order NLL. 

The NLL expression ( 2.18|) then leads to the spin-averaged form factor 

(f^ LL ) = 2(L-l-i7r)+21n(l-r 1 )(lnr 2 + z7r) + 21n(l-r 2 )(lnri + i7r) 

- ln 2 (l - n) - ln 2 (l - r 2 ) + 3 ln(l - n) + 3 ln(l - r 2 ) 

~ n) r 2 (l - r 2 ) 

l + (l-r!) 2 l + (l-r 2 ) 
with mass corrections given by ( |2.1^ ). This form is useful for comparison with other 
results on the differential cross section, as we will see in the next section. 

An analogous expression can be found for the final state emission cross section. The 
spin-averaged version of the final state form factor ( |2.36| ) is 

j" LL = 47r J BYFs(s / ,m / ) + 2(L-l-Z7r)-21nr 3 ln(l-r4) 

- 21nr 4 ln(l - r 3 ) - ln(l - r 3 ) - ln(l - r 4 ) - 2Sp(r 3 ) - 2Sp(r 4 ) 

(3.21) 



+ 2Sp(n) + 2Sp(r 2 ) + 1 ^ X + >\ l\ 2 (3.20) 



l + (l-r 3 ) 2 l + (l-r 4 ) 2 
with mass corrections given by ( |2.37| ). 
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4 Partly Differential Cross Section and Comparisons 



We may compare our spin-averaged initial state radiation form factor with one published 
in Ref. ||. In our notation, this result may be written as 

, fi N ( r i,r 2 ) + f IN (r 2 , ri ) 

J IN ~ — 7i \2 i T\ vT~ 

(1 - rx) 2 + (1 - r 2 ) z 

with 

fi N ( r i,r 2 ) = (1 -n -r 2 + rf) [4vrRe5 YFS (s, m e ) +2{L- 1)] 

^ j,„2„. , , , 2 l-n-r 2 , 2 r x 



[1 + (1 - r 2 ) 2 ] <^ \n\n) + In" — - In 



1 — r 2 1 — ri - r 2 

1 — ri — r 2 \ . 7r 2 



+ 2 Sp ( z j + 2 Sp(n + r 2 

rir 2 (l - ri - r 2 ) / 1 1 _ SrA ^ n 



1 — r 2 \ 1 — r 2 r\ J 1 — ri — r 2 

2nr 2 (l-ri-r 2 ) /(l-ri-r 2 ) 3 . 

m(l — ri — r 2 j 



ri + r 2 \ n + ^2 ri 

+ r 2 (3r 2 + 2r 1 )ln(r 1 ) - 2nr2 + + r 2 (l - r 2 ). (4.2) 

r x + r 2 1 - r 2 

The NLL limit of this expression may be obtained by summing the two collinear limits 
where r, — * separately. Carrying this out leads precisely to our spin-averaged NLL 
expression fl3.20|) . Thus, we agree with Ref. f| to NLL order. Note that the expression 
fiN does not include mass corrections, which were not calculated in Ref. ||. 

By integrating out the separate dependence on in favor of the variable z = s'/s 
= 1 — r\ — r 2 , we may obtain a result which can be compared to Ref. [f|]. To begin, we 
consider the pure real photon ISR cross section, and work in the approximation where 
the effective angles 6^- in ( |3.3| ) are replaced by a common angle 9. Then 

, ISR(O) Q2 2 

in terms of H defined by ( |3.17| ). Integrating r\ and r 2 with the constraint z — 1 — T\ — r 2 
gives 

dz 47T 2 

where r = rn 2 e (l — z)/s is the kinematic minimum value of r\ or r 2 . The result of the 
integral is, exactly, 



pi— z— ro 

cr / rfric/r 2 5(l -z-r x - r 2 )H (r 1 , r 2 ), (4.4) 



ISR(O) n2 ^ 2 

' (4-5) 



-a 



dz 47T 2 



T (l-z) 

Lj 



1 + z 2 
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Mass corrections, obtained by the prescription of Ref. [16], have the effect of replacing 
H by H + H m in §J), where 



H m (r 1 ,r 2 ) = 
Integrating the mass term gives 

"l-z-ro 



2mi ( T\ r 2 
+ ' 



r 2 T\ I 1 + z 



; H (r 1 ,r 2 ). 



(4.6) 



/ 

Jro 



dr 1 dr 2 H m (r ll r 2 ) 



1 — z s 



(4.7) 



The total mass-corrected real photon emission cross section is then 



daf R{0) qy 



dz 



4tt 5 



-(T 



(^7) {L-\) = 8*{zW 



(4i 



in the notation of Ref. |4j], where this result matches the real part of ( |2.11|) . 
The real plus virtual ISR cross section is 



da 1 



ISR(l) 



dz 



1 (Qefi 

2 V 4vr 2 



2\ 2 



O"0 



l-z-ro 



dridr 2 5(l - 2 - r x - r 2 )(f )H (r 1 ,r 2 ) 



(4.9) 



where we use the NLL expression Q3.20 ) for the virtual form factor, and will add mass 
corrections later. Doing the integral and keeping all infrared terms and terms of order L 2 
and L gives 



da 1 



ISR(l) 



dz 



4vr 2 



l-z 



L 



1+2 2 



+ K^) 2aoL (^r)i~ Lin ' +2inzin(i-z) 



+ 2Sp(l -z) + 31n2-ln 2 ^ + 



1 + z 2 



(4.10) 



where we use the notation of Ref. [HI for 



2„2 



Q2f 

4vr 2 



{27iB YFS (s,m e ) + L- 1}. 



(4.11) 



The mass correction is obtained by multiplying the pure real mass correction by the 
single virtual photon form factor evaluated in terms of s' rather than s. Thus, we add to 
the differential cross section a mass term 



da n 



2„2 



= Qle 

dr\dr 2 47r 2 



■5Y 1 {s')H m (r 1 ,r 2 ) a , 



(4.12) 
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where 

«JV) = ^{2irBYFs(s',m«)+i + lnz-l} 

- ^) + (^)ln^(21n^-L + ?-ln^). (4,3) 

Integrating over n and r 2 with z = 1 — T\ — r 2 and keeping only infrared terms and terms 
of order L, we obtain 



rfoVn _ 2,2 ( g Vtf.\ Qe e 

dz 4tt 2 1 — z 



ao{5?(s)-^Llnz}. (4.14) 



Adding the mass corrections and using the notation of ( |4.8|) gives the complete real plus 
virtual cross section at order NLL, 



tojSBQ) ...n , /^2„2\2 



dz 



+ a L{(^)ln,+ 



+ (tS) [-£lnz + 21nzln(l - *) + 31nz - In 2 z + 2 Sp(l -2)]} . (4.15) 

This result agrees precisely with the terms in (2.26) of Ref. [Q through order o?L. 
We illustrate the agreement we have found above in Figs. 2-5, for the case // = 

( 2") 

In Fig. |2|, we show the complete (3[ distribution for our exact result, our NLL and LL 
approximate results, the result of Igarashi and Nakazawa et al. 0, and the result of 
Berends et al. [|J. What we see is that there is a very good general agreement between all of 
these results. To better assess the difference between them, we plot in Fig. || the difference 
between the respective C(a 2 ) and 0(a 1 ) results. Again we see very good agreement except 
for the hardest possible photons, where then the LL result differs significantly from the 
others. 

To isolate the respective predictions for the NLL effect, we plot in Fig. |] the respective 
differences between our LL 0(a 2 ) result and the other four results. We see that there is 
again very good agreement but, at the level of 0.5 • 10 -4 , the result of Ref. [[| is somewhat 
smaller in magnitude than the other three NLL results in the Z radiative return regime 
above a cut of 0.75. 

Finally, in the Fig. |5| we isolate the size of the three NNLL results by plotting the 
difference between our NLL result and our exact result, the result from Ref. 0] and the 
result from Ref. ||J]. We again see that the most pronounced difference in the results 
occurs for the regime above v max = 0.75 where the result of Ref. jp differs by 0.5 ■ 10~ 4 
from the other two for v max < 0.975 and differs from the exact result by 1.5 • 10 -4 for 
Vmax > 0.975. The result from Ref. differs from the exact result by 0.2 • 10 -4 for 
Vmax > 0.975 but is essentially indistinguishable from it for smaller values of v max . We 
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Figure 2: This is the (3f ] distribution for the YFS3ff MC (YFS3ff is the EEX3 matrix element 
option of the /C/C MC in Ref. @), as a function of energy cut f max - It is divided by the Born 
cross-section. The IN result is from Ref. 0, and the BVNB result is from Ref. [f|. 
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Figure 3: Difference /?J 2) - for the YFS3ff MC (the EEX3 option in the KK MC), as a 
function of the cut f max - It is divided by the Born cross-section. 
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Figure 4: Next-to-leading-log contribution /?f } - (3^ L for the YFS3ff MC (the EEX3 option 
of the /C/C MC), as a function of the cut f max - It is divided by the Born cross-section. 
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Figure 5: Sub-NLL contribution (3f ] - /3^ll for the YFS3ff MC (the EEX3 option of the 
/C/C MC), as a function of the cut v max . It is divided by the Born cross-section. 
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conclude that our exact result for the 0(a 2 ) correction j3[ has a total precision tag of 
1.5 • 10~ 4 . Its NLL effect has already been implemented in the K,K MC in Ref. 0. 

We have also made the analogous study to Figs. 2-5 for 500 GeV. We find very similar 
results, with the total precision tag of 2 ■ 10~ 4 . 

5 Conclusions 

In this paper, we have presented exact results for the virtual correction to the process 
e + e _ — > ff + 7 for the ISR © FSR. The results are already in use in the JCK, MC in 
Ref. in connection with the final LEP2 data analysis. 

We have compared our results with those in Refs. |4|,|| and in general we find very 
good agreement, both at 200 GeV and at 500 GeV. For example, the size of the NNLL 
correction is shown to be at or below the level of 2 • 10~ 4 for all values of the energy cut 
parameter. Our results are fully differential and are therefore ideally suited for MC event 
generator implementation. This has been done in the /C/C MC in Ref. @. To compare our 
results with the results in Refs. @] , we have partially integrated them accordingly. While 
the results in Ref. || are also fully differential, they lack the complete mass corrections 
that our results do have. In this way, one sees that our results are in fact unique. They 
are an important part of the complete 0(a 2 ) corrections to the 2f production process 
needed for precision studies of such processes in the final LEP2 data analysis and in the 
future TESLA/LC physics. 
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Note Added 

After we had submitted this paper, we became aware of related work by G. Rodrigo, A. 
Gehrmann-De Ridder, M. Guilleaume and J. H. Kuhn, |hep-ph/0106132| . These authors 
also agree with the analogous results of Ref. |4[], when the photon azimuthal angle is 
integrated and the photon polarization is summed for the ISR process. 
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A Scalar Integrals 



Previously, in Ref. [ 10 1 , the analogous exact virtually corrected photon cross sections 



were expressed in the t channel using scalar integral functions which were calculated by 



a numerical package FF described in Ref. [|15[]. In the present paper, we have expressed 



these functions directly in terms of logarithms and dilogarithms. This appendix will give 
the expressions for the individual scalar integrals in the s channel. 

The notation for the scalar integrals will match Ref. [[l(J], with kinematic notation 
defined in Sections [l] and |^ of this paper. The scalar integrals with two denominators are 
denoted B, and the ones appearing in the form factors are 

B 12 = B(ml;m 7 ,m e ), 

B[i = B(m 2 e - sri] m 7 ,m e ), 

B 23 = B(m^;m e ,m e ), 

B 24 = B(s;m e ,m e ), 

B 34 = B(s';m e ,m e ), (A.l) 

where the first argument is the square of the momentum through the diagram, and the 
remaining arguments are the masses of the two lines. These functions are UV divergent, 
but only the following finite combinations are needed: 

/ 2 

B^-B u = In— + 6 ln^-ivr, (A.2) 

sri m A e — sri 

B 2 4~B 34 = In-, (A.3) 

s 

B 12 -B 34 = Ib-^-ztt. (A.4) 

mi 



The mass term in (|A.2| ) has been dropped when applying this expression, since the mass 



corrections are added explicitly to the massless limit of the calculation using the prescrip- 
tion of Ref. 0. 

The integrals for three and four denominators are obtained from the appendix of 
Ref. ||. For three denominators, we need the expressions 

C123 = C(mg,m^,mg- sr i ;m 1 ,m e ,m e ), 

c iu = c ( m l- sr 'i ) s / ) Wg;m 7 ,m e ,m e ) ) 

C234 = C(m^,s',s;m e ,m e ,m e ), 

C124 = C(ro^s,m^;m 7 ,me,m e ), 
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where the first three arguments are the squares of the external momenta, and the next 
three are the masses of the three lines, in cyclic order. The results are 



*rA r 23 = - ln^-Sp ll-S. --, (A.5) 

2 STi V STi J 6 



r« _ 1 1„2 S ' 1 i„2 (1 - 1 , 2 1 " T 3 



[l-r^sClL = -hi" — --ln^ - - In 

3) 134 2 m 2 m 2 r; 

s' /l — r,\ , m 2 /l — r,\ _ /m 2 
+ 2 In — In i + In — ^ In ^ - Sp — £ 

\ r i J STi \ n J \STi 

2Spf— ^-)-^-7rzln4-27rilnf^y (A.6) 



— In 2 — _. 

2 7?T.g 2 7TT,g s 

1 , n s , s , ml 2n 2 . , m 2 



{r 1 + r 2 )sC 2 u = -In 2 — --In 2 — + vriln-, (A.7) 



sC 124 = ^ln 2 — -In — ln-J- — +7riln-^. (A.8) 
2 mg m l e m z e 3 s 

The expression (|A.5|) is an analytic continuation of the one in Ref. 0, as required since 
Ti < m 2 is possible, and the expression (|A.6| ) drops mass terms. In particular, the photon 
mass regulator m 7 is dropped whenever possible. 
The expression with four denominators is 

sViDj^ = In 2 — - 2 In — In -i + In -1 In — - z'tt 

771 g ?7lg 77?,g ?7lg \ 771 g 

57T 2 s' 

+ 2Sp(r x + r 2 ) - — -2™ In — . (A.9) 

All of these expressions have been checked for agreement with the FF package. The 
function R defined in equations (|2.13| ) and ( 2.14 ) and appearing in the virtual photon 



factors is the IR-finite combination 

R(r h rj) = s (C 124 + sr 3 D r ^ u ) - srfi'^ - (1 - r t )s C[| 4 + (n + r 2 )s C 234 . (A.10) 
This completes the appendix. 
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Figure Captions 

Fig. [l| Initial state radiation graphs for e + e~ 
for / 7^ e. 



// with one virtual and one real photon, 



Fig. §. This is the (3[ > distribution for the YFS3ff MC (YFS3ff is the EEX3 matrix 
element option of the KK MC in Ref. §), as a function of energy cut v max - It is divided 
by the Born cross-section. 

Fig. |. Difference p[ 2) - (3[ 1] for the YFS3ff MC (the EEX3 option in the KK MC), as a 
function of the cut f max - It is divided by the Born cross-section. 

Fig. |. Next-to-leading-log contribution (][ 2) - f3$ h for the YFS3ff MC (the EEX3 option 
of the KK MC), as a function of the cut v max . It is divided by the Born cross-section. 

Fig. |. Sub-NLL contribution (3? - /jJJJll for the YFS3ff MC (the EEX3 option of the 
/C/C MC), as a function of the cut f max - It is divided by the Born cross-section. 
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